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Abstract 

In this paper we established the condition for a curve to satisfy 
stochastic generahzed fractional HP (Hamilton-Pontryagin) equations. 
These equations are described using Ito integral. We have also con- 
sidered the case of stochastic generalized fractional Hamiltonian equa- 
tions, for a hyperregular Lagrange function. From the stochastic ge- 
neralized fractional Hamiltonian equations, Langevin generalized frac- 
tional equations were found and numerical simulations were done. 
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1 Introduction 



J.M. Bismut was the first one who introduced concepts of stochastic ge- 
ometric mechanics, in his article from 1981, when he defined the notion of 
"stochastic Hamiltonian system". He showed that that the stochastic fiow of 
certain randomly perturbed Hamiltonian systems on flat spaces cxtrcmizes 
a stochastic action, and using this property, he proved symplecticity and the 
Noether theorem for stochastic Hamiltonian systems. Since then, there has 
been a need in to find out tools and algorithms for the study of this kind 
of systems with uncertainty. Bismut's work was continued by Lazaro-Cami 
and Ortega ([11], [12]), in the sense that his work was generalized to mani- 
folds. Stochastic Hamiltonian systems on manifolds extremize a stochastic 
action on the space of manifold valued semimartingales, the reduction of 
stochastic Hamiltonian system on the cotangent bundle of a Lie group, a 
counter example for the converse of Bismut's original theorem. Very im- 
portant in many scientific domains is fractional calculus: fractional deriva- 
tives, fractional integrals, of any real or complex order. Fractional calculus 
is used when fractional integration is needed. It is used for studying sim- 
ple dynamical systems, but it also describes complex physical systems. For 
example, applications of the fractional calculus can be found in chaotic dy- 
namics, control theory, stochastic modelling, but also in finance, hydrology, 
biophysics, physics, astrophysics, cosmology, economics and so on ([2], [4], 
[5], [9], [10]). But some other fields have just started to study problems from 
fractional point of view. It is very fashionable to study the fractional prob- 
lems of the calculus of variations and Euler-Lagrange type equations. The 
most famous fractional integrals are Riemann-Liouville, Caputo, Grunwald- 
Letnikov and the most frequently used is the Riemann-Liouville fractional 
integral. The study of Euler-Lagrange fractional equations was continued by 
Agrawal ([1], [6], [8]) that described these equations using the left, respec- 
tively right fractional derivatives in the Riemann-Liouville sense. Standard 
mult i- variable variational calculus also has some limitations. But in [13], 
C.Udriste and D. Opris showed that these limitations can by broken using 
the multi-linear control theory. In [7] the novel concepts of fractional action- 
like variational approach (FALVA) with time-dependent fractional exponent 
and exponential time-dependent term is introduced. In this paper, we re- 
strict our attention to stochastic generalized fractional Hamiltonian systems 
characterized by Wiener processes and assume that the space of admissible 
curves in configuration space is of class C^. Random effects appear in the 



2 



balance of momentum equations, as white noise, that is why we may consider 
randomly perturbed mechanical systems. It should be mentioned that the 
ideas in this paper can be readily extended to stochastic Hamiltonian sys- 
tems driven by more general semimartingales, but for the sake of clarity we 
restrict ourselves to Wiener processes. In this paper we use the generalized 
left fractional Riemann-Liouville integral defined as a mixture of the fractal 
action from physics and the discounted action at rate p, given in [7]. Within 
this context, the results of the paper are as follows: 

1. The paper presents the results from [3] which show that almost surely 
that a curve satisfies stochastic HP equations if and only if it extremizes a 
stochastic action. Suggestive examples and numerical simulations are done. 

2. Generalized fractional HP equations are described using the general- 
ized fractional Riemann-Liouville integral and the fractional Ito integral; 

3. Langevian type stochastic generalized fractional equations are obtained 
in the case of a hyperregular Lagrange function. Relevant examples and 
numerical simulations are presented. 

The paper is organized as follows: In Section 2, Hamilton-Pontryagin 
(HP) principle is given to the stochastic setting to prove that a class of me- 
chanical systems with multiplicative noise appearing as forces and torques 
possess a variational structure. For a hyperregular Lagrange function, we 
get the stochastic Hamiltonian equations that lead to Langevin equations. 
Examples and numerical simulations for the Lagrangian which describes the 
Samuelson model from economics [5] are given. In Section 3, we extend the 
generalized fractional Hamilton-Pontryagin (HP) principle to the stochastic 
setting to prove that a class of mechanical systems with multiplicative noise 
appearing as forces and torques possesses a variational structure. For a hy- 
perregular Lagrange function, we obtain the stochastic generalized fractional 
Hamiltonian equations that lead to Langevin generalized fractional equa- 
tions. For a Lagrange function, defined on iR^, the corresponding general- 
ized fractional Langevin equations are simulated. The generalized fractional 
Hamiltonian and the Lagrangian description are joined together to get the 
generahzed fractional HP system. 

2 Stochastic HP mechanics 

In this section a variational principe is introduced for a class of stochastic 
Hamilton systems on manifolds. The stochastic action is a sum of the classical 
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action and stochastic integral. The key feature of this principle is that one can 
recover stochastic Hamilton equations for these systems. Roughly speaking, 
this is accomplished by means of taking variations of this action within the 
space of curve only (not the probability space) and imposing the condition 
that this partial differential of the action must be zero. 

Let a paracompact, configuration manifold Q and J^{IR,Q) — M x TQ, 
T*Q, the associate bundle of Q. Let CP, P) be a probability space and 
{w{t), 3^t)t£[a,b], where [a, b] C M, w(t) is a real-valued Wiener process and 3^t 
is the filtration generated by the Wiener process [3]. 

The paper adopts an HP viewpoint to develop a Lagrangian description of 
stochastic Hamiltonian systems [3]. The HP principle unifies the Hamiltonian 
and Lagrangian descriptions of mechanical system. The classical HP action 
integral will be perturbed using deterministic function j : Q M. 

We consider the Lagrangian £ : J^{]R, Q) — > IR. In the stochastic context 
the HP principle states the following critical point condition on J^{]R, Q) © 
T*Q for stochastic HP action integral given by [3]: 



A{q,v,p)= I [£(s,g(s),t'(s))+<p(s),— -'y(s) >]ds+ / -i{q{s))dw{s) 



A-.nx C{PQ) M 

C{PQ) = {{s,q,v,p) e C\[a,h],PQ)\q G C\[a,hlQ),q{a) = = q,}, 

[a, h] C iR, ga, qb e Q- 



The action integral in the above principle consists of one Lebeque integral 
with respect to s and ltd stochastic integral with respect to w. The action 
is random; i.e. for every sample point a; e Q we will obtain a different time- 
dcpcndcnt Lagrangian system. We will use the following notation q(uj. s) — 
q{s), v(cl!, s) = v(s), p(uj, s) = p(s). The HP path space is a smooth infinite 
dimensional manifold. One can show that is tangent space at c = {q,v,p) G 
C{[a,b],qi,q2) consists of maps w = {q,v,p,dq,dv,dp) G C'^{[a,b],T{PQ)) 
such that 5q{a) — 5q{b) = and q,Sq are of class C^. Let {q,v,p){-,e) G 
C{PQ) denote a one-parameter family of curves in G that is differential with 
respect to e. Define the differential of A as 





(1) 



where 



dA{5q, Sv, 5p) 



d_ 
We 



■A{(jj, q{s, e),v{s, e),p{s, e))\e=o 
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where 

d 

Sq{s) = -Q^q{s,s)\e=o,Sq{a) = dq{b) = 0, 

d d 
Sv{s) = ■Q^v{s,e)\e=o,Sp{s) = —p{s,e)\e=Q. 

In terms of this differential one can state the following critical point con- 
dition: 

Theorem 1. [3] Let L : J^{]R, Q) ^ M be a Lagrangian on J^{]R, Q) of 
class with respect to t, q, v and with the globally Lipschitz first derivative 
with respect to t, q and v. Let ^ : Q ^ IR be a class function and 
with the globally Lipschitz first derivative. Then almost certainly a curve 
c— {q,v,p) e C{PQ) satisfies the stochastic HP equations: 

dq' = v'ds, 

dPi = -Q^ids+—^dw{s) 

if and only if it is a critical point of the function A : G{PQ) — > M, i.e. 

dA(c)=0. 

Let L : J^{]R, Q) ^ M he a. Lagrangian on J^{]R, Q), hyperregular, that 

/ \ 
means det I . ) ^ 0. 

From (2) the following propositions are obtained: 

Proposition 1. (Stochastic Hamilton equations). The equations 
(2) are equivalent to the following equations: 

dq = -^ds, 

where H — piV^ — L(t, q, v). 

The equations (3) represent Lagevin equations. 

Proposition 2. If L — -gijV^v^ , where gij are the components of a 
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metric on the manifold Q, equations (2) take the form: 
dq^ — v^dt, dpi — —QijV^ds H — ^^^dw{s), 



dv' = -Vji^v^v'^ds + g''^^dw{s),i,j = l..n, , 



where F*-^ are Cristoffel coefficients associated to the considered metric. 
The equations (3) become: 

dq" = g'hjds, 

Proposition 3. If L : J^{]R,1R^) M is given by: 
C, = lSijv'v^ -V{q),qeR\ 

the equations (2) take the form: 

dq"- = v^ds, 

dV d-f 

The equations (3) become: 

dq^ = S'-^pjds, 

Proposition 4. If L : J^{IR,IR'^) ^ M is given by: 

L = e-P'L{q,v),qe 
the equations (2) take the form: 

dq' = v'ds, 

Pi = e-"—. 



Proposition 5. (Samuelson. [5]) If L : J^{1R,1R^) ^ IR is given by: 
L = ~e-P'{v^ + 2avq + q^),qe M, (11) 

the equations (2) take the form: 

dq = vds, 

dp= -e'f"{av + q)ds + ^^dw{s), (12) 

p = —e~^^{v + aq). 
The equations (12) become: 

dq = —{aq + e^^p)ds, 

dp= {{a^ -l)e-f"q + ap)ds + ^. ^^^^ 

dq 

If V{q) = g^, the Euler scheme for (13) is: 
q{n + 1) = q{n) — h{aq{n) + e'^pin)) 

p(n + 1) = p{n) + h{{a^ - l)e~'^q{n) + ap{n)) + q{n)G{n),n = 0..N - 1, 

(14) 

T 

where T > 0, h = — , > 0, G{n) = w{{n + l)h) — w{nh) and q{n) = 

q{uj, nh), p{n) = p{uj, nh), p E (0, 1), a E (—1, 1). 

For p = 0.003, a = 0.03, h = 0.001 using Maple 13, the orbit {n,q{nh)) 
is represented in Fig 1 and (n, q{uj, nh)) in Fig 2: 
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In Figures 3 and 4 we can visualize the orbits {n,p{nh)), {n,p{uj,nh)): 



Fig 3. {n,p{nh)) Fig 4. {n,p{uj,nh)) 




Figures 5 and 6 represent the orbits {q{nh),p{nh)) and (g(a;, nh),p{u!, nh)): 



Fig 5. {q{nh) , p{nh)) Fig 6. {q{uj , nh) , p{uj , nh)) 




3 Stochastic generalized fractional HP 

In this section a generalized variational principle is introduced for a class 
of stochastic generalized fractional Hamiltonian systems on manifold. We 
use the generalized left fractional Riemann-Liouville integral ([7]) defined as 
a mixture of the fractal action from physics and the discounted action at rate 
p. Let / : iR — > iR be an integrable function, a : M ^ M a. function of class 
C^. The generalized left fractional Riemman-Liouville integral is given by: 
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= f ^ , / ^.J {s){t - 5)"(-*)-^e^(-*)rf5 (15) 



where 



poo 

Vi{a{s-t)) =V{a{z))\z^ s-t,T{a{z)) = / (s - t)"(")-^e-("-*)(is, 

JQ 

is the modified Euler Gamma function. 

If a{z) = a = const., < a < 1, p = 0, from (15) we obtain the left 
fractional Riemman-Liouville integral ([6], [7], [8]). In fact, the generalized 
left fractional Riemann-Liouville integral is a generalization of the single time 
Stieltyes integral ([13]). 

In (15), s is the intrinsec time and t is the observer time, t ^ s. Let 
g : IR^ Mhe the function: 

ri{a[s - 1)) 

We consider £ : {M, Q) ^ M and j : Q ^ JR. In the stochastic context 
the HP principle states the following critical point condition on J^{1R, Q) © 
T*{Q) for stochastic HP generahzed fractional action given by: 

, 

A''{q,v,p,t)= / [£(s,g(s),v(s))+<p(s),-^-v(s) >]gt{s)d»^ 

Jto ds ^^^^ 



+ / l{l{s))gt{s)dw{s). 



to 

The first integral in (17) is a Lebeques integral with respect to s and the 
second one is an Ito integral. 
Using Theorem 1, we get: 

Theorem 2. If L : J^{1R, Q) ^ M and : Q ^ M satisfy the hypothesis 
from Theorem 1, then almost certainly a curve c = {q,v,p) e C{PQ) satisfies 
the stochastic HP equations with intrinsec and observer times: 
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dg* = v'ds, 

,da{s — t), , , a(s — t) — l, 

1 



Pi = = / s. 

From (18) wc obtain: 

(i) If a{z) = 1 and p = then equations (2) are obtained; 

(ii) If a{z) = a = const., 0<a<l, p = then the following relations 
are deduced from (18): 



dq^ ^ v^ds, 

= W "^'7^^^' + (19) 
Pi = = l..n,t / s. 

The equations (19) represent the stochastic fractional equations. 
If L : J^{]R,Q) — > iR is hyperregular, using (18) the following proposi- 
tions hold: 

Proposition 6. (Stochastic generalized fractional Hamilton equations.) 
Equations (18) are equivalent with the equations: 



dq = -^ds, 

OPi 

dpi = ~ PiK^^ t))ds + ^^Q^} dw{s), 



(20) 



where 



H ^Piv' - L{s,q,v), 

dais — t) , , , a(s — t) — 1 1 dTiiais — t)) 

h(s, t) = — ^ -Int -s + ^ hp- 7-^^ TT — -■ 

^ ' ^ ds ' ' s-t Vi{a{s-t)) ds 

(21) 
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The equations (20) represent the generahzed fractional Langevin equations. 

Proposition 7. If L — -Qijv'^v-' , where gij are the components of a 
metric on the manifold Q, then equation (18) takes the form: 



dq^ ^ v^ds, 

■ u -diiq) (22) 

dv^ = -{r.^^v^v^ - h{s,ty)ds + g^^^dw{s),i,j = l..n, 

where V^/^ are Cristoffel coefficients associated to the considered metric and 
h{s,t) is given by (21). The equations (20) become: 

dq^ = g^^pjds, 

dp. = il^P^p' - his,t)p,)ds + ^dwisU^ l..n, ^''^ 

where h{s,t) is given by (21). 

The equations (23) can be used for generahzed fractional motion of rela- 
tivistic particle with noise ([7]). 

Proposition 8. // £ : {M, M) ^ Ft is given by: 

i:{q,v)^^-v'-V{q) (24) 
where V : M ^ M and j : M ^ M, then equations (23) are given by: 



dq — pds, 

dV d-f(q) (25) 

dp — ( — ; h(s, t)p)ds H ; — dwis). 

dq dq 

If V{q) — cos{q), 7(g) = sin{q), the Euler scheme for (25) is given by: 

q{n + 1) = q{n) + kp{n)) 

p{n + 1) = p{n) + k{sin{q{n)) — h{nk, t j) + cos{q{nj)G{n), n = 0..N — 1, 
T 

where T > 0, k = —, N > 0, G{n) = w{{n + l)h) - w{nh) and 
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, , , , 1 / , N, I ,1 a(nk — t) — 1 

n[nk, t) = a ink, tjlnlt — nk\ H ; h p— 

nk — t 

r\a{nk~t)), 



Ti{a{nk — t)) 

a\s -t) = T\a{s - t)) = d^^i^if-')), 

as as 

For a{s — t) = a, a = 0.6, k = 0.001, t = 0.8 using Maple 13 the orbit 

{n,p{nk)) is represented in Figure 7 and the orbit {n,p{u, nk)) is represented 

in Figure 8: 



Fig 7. {n,p{nk)) Fig 8. {n,p{u,nk)) 




Figure 9 displays the orbit {q{nk) , p{nk)) and Figure 10 displays the orbit 

{q{uj, nk),p{u, nk)) 



Fig 9. {q{nk) , p{nk)) Fig 10. {q^uj , nk) , p{uj , nk)) 
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4 Conclusions 



In this paper we have described the stochastic generahzed fractional HP 
principle, using the classical stochastic HP principle [3]. Using a hyperregular 
Lagrange function, Langevin-type generalized fractional equations were illus- 
trated. Wc have done the numerical simulations for the case of a Lagrange 
function defined on IR^. In our future papers we will study the stochastic 
stability of the obtained equations al well as the description of the credibility 
generalized fractional HP principle. 
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